FGCD: a database of fungal gene clusters related to secondary metabolism.

Fungal secondary metabolites are not necessary for growth, but they are important for fungal metabolism and ecology because they provide selective advantages for competition, survival and interactions with the environment. These various metabolites are widely used as medicinal precursors and insecticides. Secondary metabolism genes are commonly arranged in clusters along chromosomes, which allow for the coordinate control of complete pathways. In this study, we created the Fungal Gene Cluster Database to store, retrieve, and visualize secondary metabolite gene cluster information across fungal species. The database was created by merging data from RNA sequencing, Basic Local Alignment Search Tool, genome browser, enrichment analysis and the R Shiny web framework to visualize and query putative gene clusters. This database facilitated the rapid and thorough examination of significant gene clusters across fungal species by detecting, defining and graphically displaying the architecture, organization and expression patterns of secondary metabolite gene clusters. In general, this genomic resource makes use of the tremendous chemical variety of the products of these ecologically and biotechnologically significant gene clusters to our further understanding of fungal secondary metabolism. Database URL: https://www.hebaubioinformatics.cn/FungalGeneCluster/.

Genomic insights into Penicillium chrysogenum adaptation to subseafloor sedimentary environments.

BACKGROUND: Penicillium chrysogenum is a filamentous fungal species with diverse habitats, yet little is known about its genetics in adapting to extreme subseafloor sedimental environments. RESULTS: Here, we report the discovery of P. chrysogenum strain 28R-6-F01, isolated from deep coal-bearing sediments 2306 m beneath the seafloor. This strain possesses exceptional characteristics, including the ability to thrive in extreme conditions such as high temperature (45 °C), high pressure (35 Mpa), and anaerobic environments, and exhibits broad-spectrum antimicrobial activity, producing the antibiotic penicillin at a concentration of 358 µg/mL. Genome sequencing and assembly revealed a genome size of 33.19 Mb with a GC content of 48.84%, containing 6959 coding genes. Comparative analysis with eight terrestrial strains identified 88 unique genes primarily associated with penicillin and aflatoxins biosynthesis, carbohydrate degradation, viral resistance, and three secondary metabolism gene clusters. Furthermore, significant expansions in gene families related to DNA repair were observed, likely linked to the strain's adaptation to its environmental niche. CONCLUSIONS: Our findings provide insights into the genomic and biological characteristics of P. chrysogenum adaptation to extreme anaerobic subseafloor sedimentary environments, such as high temperature and pressure.

Analysis of five near-complete genome assemblies of the tomato pathogen Cladosporium fulvum uncovers additional accessory chromosomes and structural variations induced by transposable elements effecting the loss of avirulence genes.

BACKGROUND: Fungal plant pathogens have dynamic genomes that allow them to rapidly adapt to adverse conditions and overcome host resistance. One way by which this dynamic genome plasticity is expressed is through effector gene loss, which enables plant pathogens to overcome recognition by cognate resistance genes in the host. However, the exact nature of these loses remains elusive in many fungi. This includes the tomato pathogen Cladosporium fulvum, which is the first fungal plant pathogen from which avirulence (Avr) genes were ever cloned and in which loss of Avr genes is often reported as a means of overcoming recognition by cognate tomato Cf resistance genes. A recent near-complete reference genome assembly of C. fulvum isolate Race 5 revealed a compartmentalized genome architecture and the presence of an accessory chromosome, thereby creating a basis for studying genome plasticity in fungal plant pathogens and its impact on avirulence genes. RESULTS: Here, we obtained near-complete genome assemblies of four additional C. fulvum isolates. The genome assemblies had similar sizes (66.96 to 67.78 Mb), number of predicted genes (14,895 to 14,981), and estimated completeness (98.8 to 98.9%). Comparative analysis that included the genome of isolate Race 5 revealed high levels of synteny and colinearity, which extended to the density and distribution of repetitive elements and of repeat-induced point (RIP) mutations across homologous chromosomes. Nonetheless, structural variations, likely mediated by transposable elements and effecting the deletion of the avirulence genes Avr4E, Avr5, and Avr9, were also identified. The isolates further shared a core set of 13 chromosomes, but two accessory chromosomes were identified as well. Accessory chromosomes were significantly smaller in size, and one carried pseudogenized copies of two effector genes. Whole-genome alignments further revealed genomic islands of near-zero nucleotide diversity interspersed with islands of high nucleotide diversity that co-localized with repeat-rich regions. These regions were likely generated by RIP, which generally asymmetrically affected the genome of C. fulvum. CONCLUSIONS: Our results reveal new evolutionary aspects of the C. fulvum genome and provide new insights on the importance of genomic structural variations in overcoming host resistance in fungal plant pathogens.

Whole genome sequencing and annotation of Aspergillus flavus JAM-JKB-B HA-GG20.

Groundnuts are mostly contaminated with the mold Aspergillus flavus which produces a carcinogenic mycotoxin called as aflatoxin. It is very important to understand the genetic factors underlying its pathogenicity, regulation, and biosynthesis of secondary metabolites and animal toxicities, but it still lacks useful information due to certain gaps in the era of modern technology. Therefore, the present study was considered to determine the key genes and metabolites involved in the biosynthesis of aflatoxin by using a molecular approach in a virulent strain of Aspergillus. The whole genome sequence of highly toxic and virulent Aspergillus isolates JAM-JKB-B HA-GG20 revealed 3,73,54,834 bp genome size, 2, 26, 257 number of contigs with N50 value of 49,272 bp, 12,400 genes and 48.1% of GC contained respectively. The genome sequence was compared with other known aflatoxin producing and non-producing genome of Aspergillus spp. and 61 secondary metabolite (SM) gene clusters were annotated with the toxic strain JAM-JKB-BHA-GG20 which showed similarity with other Aspergillus spp. A total number of eight genes (ver-1, AflR, pksA, uvm8, omt1, nor-1, Vha and aflP) were identified related to biosynthesis of aflatoxin and ochratoxin. Also, 69 SSR with forward and reverse primers and 137 di and tri nucleotide motifs were identified in the nucleotide sequence region related to aflatoxin gene pathway. The genes and putative metabolites identified in this study are potentially involved in host invasion and pathogenicity. As such, the genomic information obtained in this study is helpful in understanding aflatoxin gene producing pathway in comparison to other Aspergillus spp. and predicted presence of other secondary metabolites clusters viz. Nrps, T1pks etc. genes associated with a biosynthesis of OTA mycotoxin.

A rapid PCR-based method to determine the Neurospora crassa mating type.

So far mating type determination in Neurospora crassa requires test crosses with strains of known mating type. We present a simple, quick, and reliable polymerase chain reaction-based method for mating type determination in N. crassa.

Multifunctionality of AsCFEM6 and AsCFEM12 effectors from the potato early blight pathogen Alternaria solani.

Plant pathogens secrete fungal-specific common in several fungal extracellular membrane (CFEM) effectors to manipulate host immunity and contribute to their virulence. Little is known about effectors and their functions in Alternaria solani, the necrotrophic fungal pathogen causing potato early blight. To identify candidate CFEM effector genes, we mined A. solani genome databases. This led to the identification of 12 genes encoding CFEM proteins (termed AsCFEM1-AsCFEM12) and 6 of them were confirmed to be putative secreted effectors. In planta expression revealed that AsCFEM6 and AsCFEM12 have elicitor function that triggers plant defense response including cell death in different botanical families. Targeted gene disruption of AsCFEM6 and AsCFEM12 resulted in a change in spore development, significant reduction of virulence on potato and eggplant susceptible cultivars, increased resistance to fungicide stress, variation in iron acquisition and utilization, and the involvement in 1,8-dihydroxynaphthalene (DHN) melanin biosynthesis pathway. Using maximum likelihood method, we found that positive selection likely caused the polymorphism within AsCFEM6 and AsCFEM12 homologs in different Alternaria spp. Site-directed mutagenesis analysis indicated that positive selection sites within their CFEM domains are required for cell death induction in Nicotiana benthamiana and are critical for response to abiotic stress in yeast. These results demonstrate that AsCFEM effectors possess additional functions beyond their roles in host plant immune response and pathogen virulence.

Complexity of fungal polyketide biosynthesis and function.

Where does one draw the line between primary and secondary metabolism? The answer depends on the perspective. Microbial secondary metabolites (SMs) were at first believed not to be very important for the producers because they are dispensable for growth under laboratory conditions. However, such compounds become important in natural niches of the organisms, and some are of prime importance for humanity. Polyketides are an important group of SMs with aflatoxin as a well-known and well-characterized example. In Aspergillus spp., all 34 afl genes encoding the enzymes for aflatoxin biosynthesis are located in close vicinity on chromosome III in a so-called gene cluster. This led to the assumption that most genes required for polyketide biosynthesis are organized in gene clusters. Recent research, however, revealed an enormous complexity of the biosynthesis of different polyketides, ranging from individual polyketide synthases to a gene cluster producing several compounds, or to several clusters with additional genes scattered in the genome for the production of one compound. Research of the last decade furthermore revealed a huge potential for SM biosynthesis hidden in fungal genomes, and methods were developed to wake up such sleeping genes. The analysis of organismic interactions starts to reveal some of the ecological functions of polyketides for the producing fungi.

Building a eukaryotic chromosome arm by de novo design and synthesis.

The genome of an organism is inherited from its ancestor and continues to evolve over time, however, the extent to which the current version could be altered remains unknown. To probe the genome plasticity of Saccharomyces cerevisiae, here we replace the native left arm of chromosome XII (chrXIIL) with a linear artificial chromosome harboring small sets of reconstructed genes. We find that as few as 12 genes are sufficient for cell viability, whereas 25 genes are required to recover the partial fitness defects observed in the 12-gene strain. Next, we demonstrate that these genes can be reconstructed individually using synthetic regulatory sequences and recoded open-reading frames with a "one-amino-acid-one-codon" strategy to remain functional. Finally, a synthetic neochromsome with the reconstructed genes is assembled which could substitute chrXIIL for viability. Together, our work not only highlights the high plasticity of yeast genome, but also illustrates the possibility of making functional eukaryotic chromosomes from entirely artificial sequences.

Lineage-specific genes are clustered with HET-domain genes and respond to environmental and genetic manipulations regulating reproduction in Neurospora.

Lineage-specific genes (LSGs) have long been postulated to play roles in the establishment of genetic barriers to intercrossing and speciation. In the genome of Neurospora crassa, most of the 670 Neurospora LSGs that are aggregated adjacent to the telomeres are clustered with 61% of the HET-domain genes, some of which regulate self-recognition and define vegetative incompatibility groups. In contrast, the LSG-encoding proteins possess few to no domains that would help to identify potential functional roles. Possible functional roles of LSGs were further assessed by performing transcriptomic profiling in genetic mutants and in response to environmental alterations, as well as examining gene knockouts for phenotypes. Among the 342 LSGs that are dynamically expressed during both asexual and sexual phases, 64% were detectable on unusual carbon sources such as furfural, a wildfire-produced chemical that is a strong inducer of sexual development, and the structurally-related furan 5-hydroxymethyl furfural (HMF). Expression of a significant portion of the LSGs was sensitive to light and temperature, factors that also regulate the switch from asexual to sexual reproduction. Furthermore, expression of the LSGs was significantly affected in the knockouts of adv-1 and pp-1 that regulate hyphal communication, and expression of more than one quarter of the LSGs was affected by perturbation of the mating locus. These observations encouraged further investigation of the roles of clustered lineage-specific and HET-domain genes in ecology and reproduction regulation in Neurospora, especially the regulation of the switch from the asexual growth to sexual reproduction, in response to dramatic environmental conditions changes.

A fungal sesquiterpene biosynthesis gene cluster critical for mutualist-pathogen transition in Colletotrichum tofieldiae.

Plant-associated fungi show diverse lifestyles from pathogenic to mutualistic to the host; however, the principles and mechanisms through which they shift the lifestyles require elucidation. The root fungus Colletotrichum tofieldiae (Ct) promotes Arabidopsis thaliana growth under phosphate limiting conditions. Here we describe a Ct strain, designated Ct3, that severely inhibits plant growth. Ct3 pathogenesis occurs through activation of host abscisic acid pathways via a fungal secondary metabolism gene cluster related to the biosynthesis of sesquiterpene metabolites, including botrydial. Cluster activation during root infection suppresses host nutrient uptake-related genes and changes mineral contents, suggesting a role in manipulating host nutrition state. Conversely, disruption or environmental suppression of the cluster renders Ct3 beneficial for plant growth, in a manner dependent on host phosphate starvation response regulators. Our findings indicate that a fungal metabolism cluster provides a means by which infectious fungi modulate lifestyles along the parasitic-mutualistic continuum in fluctuating environments.

Extending the reach of homology by using successive computational filters to find yeast pheromone genes.

The mating of fungi depends on pheromones that mediate communication between two mating types. Most species use short peptides as pheromones, which are either unmodified (e.g., α-factor in Saccharomyces cerevisiae) or C-terminally farnesylated (e.g., a-factor in S. cerevisiae). Peptide pheromones have been found by genetics or biochemistry in a small number of fungi, but their short sequences and modest conservation make it impossible to detect homologous sequences in most species. To overcome this problem, we used a four-step computational pipeline to identify candidate a-factor genes in sequenced genomes of the Saccharomycotina, the fungal clade that contains most of the yeasts: we require that candidate genes have a C-terminal prenylation motif, are shorter than 100 amino acids long, and contain a proteolytic-processing motif upstream of the potential mature pheromone sequence and that closely related species contain highly conserved homologs of the potential mature pheromone sequence. Additional manual curation exploits the observation that many species carry more than one a-factor gene, encoding identical or nearly identical pheromones. From 332 Saccharomycotina genomes, we identified strong candidate pheromone genes in 241 genomes, covering 13 clades that are each separated from each other by at least 100 million years, the time required for evolution to remove detectable sequence homology among small pheromone genes. For one small clade, the Yarrowia, we demonstrated that our algorithm found the a-factor genes: deleting all four related genes in the a-mating type of Yarrowia lipolytica prevents mating.

The Bcvic1 and Bcvic2 vegetative incompatibility genes in Botrytis cinerea encode proteins with domain architectures involved in allorecognition in other filamentous fungi.

Vegetative incompatibility is a fungal allorecognition system characterised by the inability of genetically distinct conspecific fungal strains to form a viable heterokaryon and is controlled by multiple polymorphic loci termed vic (vegetative incompatibility) or het (heterokaryon incompatibility). We have genetically identified and characterised the first vic locus in the economically important, plant-pathogenic, necrotrophic fungus Botrytis cinerea. A bulked segregant approach coupled with whole genome Illumina sequencing of near-isogenic lines of B. cinerea was used to map a vic locus to a 60-kb region of the genome. Within that locus, we identified two adjacent, highly polymorphic open reading frames, Bcvic1 and Bcvic2, which encode predicted proteins that contain domain architectures implicated in vegetative incompatibility in other filamentous fungi. Bcvic1 encodes a predicted protein containing a putative serine esterase domain, a NACHT family of NTPases domain, and several Ankyrin repeats. Bcvic2 encodes a putative syntaxin protein containing a SNARE domain; such proteins typically function in vesicular transport. Deletion of Bcvic1 and Bcvic2 individually had no effect on vegetative incompatibility. However, deletion of the region containing both Bcvic1 and Bcvic2 resulted in mutant lines that were severely restricted in growth and showed loss of vegetative incompatibility. Complementation of these mutants by ectopic expression restored the growth and vegetative incompatibility phenotype, indicating that Bcvic1 and Bcvic2 are controlling vegetative incompatibility at this vic locus.

Microbe-induced gene silencing boosts crop protection against soil-borne fungal pathogens.

Small RNA (sRNA)-mediated trans-kingdom RNA interference (RNAi) between host and pathogen has been demonstrated and utilized. However, interspecies RNAi in rhizospheric microorganisms remains elusive. In this study, we developed a microbe-induced gene silencing (MIGS) technology by using a rhizospheric beneficial fungus, Trichoderma harzianum, to exploit an RNAi engineering microbe and two soil-borne pathogenic fungi, Verticillium dahliae and Fusarium oxysporum, as RNAi recipients. We first detected the feasibility of MIGS in inducing GFP silencing in V. dahliae. Then by targeting a fungal essential gene, we further demonstrated the effectiveness of MIGS in inhibiting fungal growth and protecting dicotyledon cotton and monocotyledon rice plants against V. dahliae and F. oxysporum. We also showed steerable MIGS specificity based on a selected target sequence. Our data verify interspecies RNAi in rhizospheric fungi and the potential application of MIGS in crop protection. In addition, the in situ propagation of a rhizospheric beneficial microbe would be optimal in ensuring the stability and sustainability of sRNAs, avoiding the use of nanomaterials to carry chemically synthetic sRNAs. Our finding reveals that exploiting MIGS-based biofungicides would offer straightforward design and implementation, without the need of host genetic modification, in crop protection against phytopathogens.

Unraveling the role of effector proteins in Bipolaris oryzae infecting North East Indian rice cultivars through time-course transcriptomics analysis.

Bipolaris oryzae, causing brown spot disease in rice, is one of the neglected diseases reducing rice productivity. Limited knowledge is available on the genetics of host-pathogen interaction. Here, we used time-course transcriptome sequencing to elucidate the differential transcriptional responses of the pathogen genes in two contradictory infection-responsive rice hosts. Evaluation of transcriptome data showed similar regulation of fungal genes within susceptible (1733) and resistant (1846) hosts at an early stage however, in the later stage, the number was significantly higher in susceptible (2877) compared to resistant (1955) hosts. GO enrichment terms for upregulated genes showed a similar pattern in both the hosts at an early stage, but in the later stage terms related to degradation of carbohydrates, carbohydrate transport, and pathogenesis are enriched extensively within the susceptible host. Likewise, similar expression responses were observed with the secretory and effector proteins. Plant pathogenic homologs genes such as those involved in appressorium and conidia formation, host cell wall degradative enzymes, etc. were reported to be highly upregulated within the susceptible host. This study predicts the successful establishment of B. oryzae BO1 in both the host surfaces at an early stage, while disease progression only occurs in the susceptible host in later stage.

Application of CRISPR in Filamentous Fungi and Macrofungi: From Component Function to Development Potentiality.

Fungi, particularly filamentous fungi and macrofungi, have a very powerful ability to produce secondary metabolites and can serve as excellent chassis cells for the production of enzymes or natural products of great value in synthetic biology. Thus, it is imperative to establish simple, reliable, and efficient techniques for their genetic modification. However, the heterokaryosis of some fungi and the dominance of nonhomologous end-joining (NHEJ) repair mechanisms in vivo have been greatly affecting the efficiency of fungal gene editing. In recent years, the CRISPR/Cas9 system has been applied as a widely used gene editing technology in life science research and has also played an important role in the genetic modification of filamentous and macrofungi. The various functional components (cas9, sgRNA, promoter, and screening marker) of the CRISPR/Cas9 system and its development, as well as the difficulties and potential of the CRISPR/Cas9 system in filamentous fungus and macrofungi, are the main topics of this paper.

Comparison of the impact of two key fungal signalling pathways on Zymoseptoria tritici infection reveals divergent contribution to invasive growth through distinct regulation of infection-associated genes.

The lifecycle of Zymoseptoria tritici requires a carefully regulated asymptomatic phase within the wheat leaf following penetration of the mesophyll via stomata. Here we compare the roles in this process of two key fungal signalling pathways, mutants of which were identified through forward genetics due to their avirulence on wheat. Whole-genome resequencing of avirulent Z. tritici T-DNA transformants identified disruptive mutations in ZtBCK1 from the kinase cascade of the cell wall integrity (CWI) pathway, and the adenylate cyclase gene ZtCYR1. Targeted deletion of these genes abolished the pathogenicity of the fungus and led to similar in vitro phenotypes to those associated with disruption of putative downstream kinases, both supporting previous studies and confirming the importance of these pathways in virulence. RNA sequencing was used to investigate the effect of ZtBCK1 and ZtCYR1 deletion on gene expression in both the pathogen and host during infection. ZtBCK1 was found to be required for the adaptation to the host environment, controlling expression of infection-associated secreted proteins, including known virulence factors. Meanwhile, ZtCYR1 is implicated in controlling the switch to necrotrophy, regulating expression of effectors associated with this transition. This represents the first study to compare the influence of CWI and cAMP signalling on in planta transcription of a fungal plant pathogen, providing insights into their differential regulation of candidate effectors during invasive growth.

Permissiveness and competition within and between Neurospora crassa syncytia.

A multinucleate syncytium is a common growth form in filamentous fungi. Comprehensive functions of the syncytial state remain unknown, but it likely allows for a wide range of adaptations to enable filamentous fungi to coordinate growth, reproduction, responses to the environment, and to distribute nuclear and cytoplasmic elements across a colony. Indeed, the underlying mechanistic details of how syncytia regulate cellular and molecular processes spatiotemporally across a colony are largely unexplored. Here, we implemented a strategy to analyze the relative fitness of different nuclear populations in syncytia of Neurospora crassa, including nuclei with loss-of-function mutations in essential genes, based on production of multinucleate asexual spores using flow cytometry of pairings between strains with differentially fluorescently tagged nuclear histones. The distribution of homokaryotic and heterokaryotic asexual spores in pairings was assessed between different auxotrophic and morphological mutants, as well as with strains that were defective in somatic cell fusion or were heterokaryon incompatible. Mutant nuclei were compartmentalized into both homokaryotic and heterokaryotic asexual spores, a type of bet hedging for maintenance and evolution of mutational events, despite disadvantages to the syncytium. However, in pairings between strains that were blocked in somatic cell fusion or were heterokaryon incompatible, we observed a "winner-takes-all" phenotype, where asexual spores originating from paired strains were predominantly one genotype. These data indicate that syncytial fungal cells are permissive and tolerate a wide array of nuclear functionality, but that cells/colonies that are unable to cooperate via syncytia formation actively compete for resources.

Divergence of alternative sugar preferences through modulation of the expression and activity of the Gal3 sensor in yeast.

Optimized nutrient utilization is crucial for the progression of microorganisms in competing communities. Here we investigate how different budding yeast species and ecological isolates have established divergent preferences for two alternative sugar substrates: Glucose, which is fermented preferentially by yeast, and galactose, which is alternatively used upon induction of the relevant GAL metabolic genes. We quantified the dose-dependent induction of the GAL1 gene encoding the central galactokinase enzyme and found that a very large diversification exists between different yeast ecotypes and species. The sensitivity of GAL1 induction correlates with the growth performance of the respective yeasts with the alternative sugar. We further define some of the mechanisms, which have established different glucose/galactose consumption strategies in representative yeast strains by modulating the activity of the Gal3 inducer. (1) Optimal galactose consumers, such as Saccharomyces uvarum, contain a hyperactive GAL3 promoter, sustaining highly sensitive GAL1 expression, which is not further improved upon repetitive galactose encounters. (2) Desensitized galactose consumers, such as S. cerevisiae Y12, contain a less sensitive Gal3 sensor, causing a shift of the galactose response towards higher sugar concentrations even in galactose experienced cells. (3) Galactose insensitive sugar consumers, such as S. cerevisiae DBVPG6044, contain an interrupted GAL3 gene, causing extremely reluctant galactose consumption, which is, however, improved upon repeated galactose availability. In summary, different yeast strains and natural isolates have evolved galactose utilization strategies, which cover the whole range of possible sensitivities by modulating the expression and/or activity of the inducible galactose sensor Gal3.

Fission yeast Swi2 designates cell-type specific donor and stimulates Rad51-driven strand exchange for mating-type switching.

A haploid of the fission yeast Schizosaccharomyces pombe expresses either the P or M mating-type, determined by the active, euchromatic, mat1 cassette. Mating-type is switched by Rad51-driven gene conversion of mat1 using a heterochromatic donor cassette, mat2-P or mat3-M. The Swi2-Swi5 complex, a mating-type switching factor, is central to this process by designating a preferred donor in a cell-type-specific manner. Swi2-Swi5 selectively enables one of two cis-acting recombination enhancers, SRE2 adjacent to mat2-P or SRE3 adjacent to mat3-M. Here, we identified two functionally important motifs in Swi2, a Swi6 (HP1 homolog)-binding site and two DNA-binding AT-hooks. Genetic analysis demonstrated that the AT-hooks were required for Swi2 localization at SRE3 to select the mat3-M donor in P cells, while the Swi6-binding site was required for Swi2 localization at SRE2 to select mat2-P in M cells. In addition, the Swi2-Swi5 complex promoted Rad51-driven strand exchange in vitro. Taken together, our results show how the Swi2-Swi5 complex would localize to recombination enhancers through a cell-type specific binding mechanism and stimulate Rad51-driven gene conversion at the localization site.

Combined pangenomics and transcriptomics reveals core and redundant virulence processes in a rapidly evolving fungal plant pathogen.

BACKGROUND: Studying genomic variation in rapidly evolving pathogens potentially enables identification of genes supporting their "core biology", being present, functional and expressed by all strains or "flexible biology", varying between strains. Genes supporting flexible biology may be considered to be "accessory", whilst the "core" gene set is likely to be important for common features of a pathogen species biology, including virulence on all host genotypes. The wheat-pathogenic fungus Zymoseptoria tritici represents one of the most rapidly evolving threats to global food security and was the focus of this study. RESULTS: We constructed a pangenome of 18 European field isolates, with 12 also subjected to RNAseq transcription profiling during infection. Combining this data, we predicted a "core" gene set comprising 9807 sequences which were (1) present in all isolates, (2) lacking inactivating polymorphisms and (3) expressed by all isolates. A large accessory genome, consisting of 45% of the total genes, was also defined. We classified genetic and genomic polymorphism at both chromosomal and individual gene scales. Proteins required for essential functions including virulence had lower-than average sequence variability amongst core genes. Both core and accessory genomes encoded many small, secreted candidate effector proteins that likely interact with plant immunity. Viral vector-mediated transient in planta overexpression of 88 candidates failed to identify any which induced leaf necrosis characteristic of disease. However, functional complementation of a non-pathogenic deletion mutant lacking five core genes demonstrated that full virulence was restored by re-introduction of the single gene exhibiting least sequence polymorphism and highest expression. CONCLUSIONS: These data support the combined use of pangenomics and transcriptomics for defining genes which represent core, and potentially exploitable, weaknesses in rapidly evolving pathogens.